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We have explored the manifold physical scenario emerging from a model of Dynamic Nuclear 
Polarization (DNP) via thermal mixing under the hypothesis of highly effective electron-electron 
interaction. When the electron and nuclear reservoirs are also assumed to be in strong thermal 
contact and the microwave irradiation saturates the target electron transition, the enhancement of 
the nuclear polarization is expected to be considerably high even if the irradiation frequency is set far 
away from the centre of the ESR line (as already observed by Borghini) and the typical polarization 
time is reduced on moving towards the boundaries of said line. More reasonable behaviours are 
obtained by reducing the level of microwave saturation or the contact between electrons and nuclei 
in presence of nuclear leakage. In both cases the function describing the dependency of the steady 
state nuclear polarization on the frequency of irradiation becomes sharper at the edges and the build 
up rate decreases on moving off-resonance. If qualitatively similar in terms of the effects produced 
on nuclear polarization, the degree of microwave saturation and of electron-nucleus contact has 
a totally different impact on electron polarization, which is of course strongly correlated to the 
effectiveness of saturation and almost insensitive, at the steady state, to the magnitude of the 
interactions between the two spin reservoirs. The likelihood of the different scenario is discussed in 
the light of the experimental data currently available in literature, to point out which aspects are 
suitably accounted and which are not by the declinations of thermal mixing DNP considered here. 



I. INTRODUCTION 

In the last decade Dynamic Nuclear Polarization 
(DNP) has established itself as a powerful technique to 
overcome the limited sensitivity of Nuclear Magnetic Res- 
onance (NMR) [1]. More recently, as a consequence 
of the impressive experimental results over a wide area 
of applications, ranging from analytical [H [3] to poten- 
tially diagnostic methods [IHS], the scientific community 
has started to deepen the existing theoretical knowledge 
about the physics of the polarization process [THS]. De- 
pending on the specific conditions of the experiment, the 
transfer of magnetic order from the electron to the nu- 
clear system occurs by different mechanisms, named Solid 
Effect dnmi], Cross Effect (HHTl] and Thermal Mixing 
[TTl [T51 [T7] . The latter regime is believed to apply to 
those samples and experimental conditions typically ex- 
ploited in biomedical applications [18], which nowadays 
are attracting a large interest. 

The original theoretical description of low temperature 
DNP via Thermal Mixing (TM) is due to Borghini [T6^ 
and re-proposed in a slighly different fashion by Abragam 
and Goldman in their famous review The model is 
based on the hypothesis oi (i) very efficient spectral dif- 
fusion, (ii) complete saturation of the irradiated Electron 
Spin Resonance (ESR) isocromate and (Hi) existence of 
a perfect contact between electrons and nuclei. This lat- 
ter forcing the establishment of a common temperature 
between nuclear and electron reservoirs at any time. 

Despite the Borghini prediction qualitatively depicts 
some aspects of the experimental scenario, no informa- 
tion about the dynamics of the process is given, while the 
steady state nuclear polarization is always overestimated. 



The quantitative agreement is especially poor when mov- 
ing from the centre to the edges of the ESR spectrum. 
In order to reduce the discrepancies between theory and 
experimental observations in the nuclear steady state be- 
haviour, Jannin et al. |19j has recently proposed a vari- 
ant of the Borghini model where the irradiated ESR line 
portion is only partially saturated. Again, the dynamical 
problem has not been tackled. 

A general methodology to compute the full time evo- 
lution of the nuclear polarization in the low temperature 
TM regime, relying on a mean field approach and based 
on a proper system of rate equations, has been described 
in |20j. In this work we exploited that mathematical 
treatment (briefly recalled in Section |ll]) for providing a 
comprehensive picture of the electron and nuclear polar- 
ization dynamics (including the relevant steady states), 
over the whole microwave spectrum, for different choices 
of the five time constants describing the basic interactions 
and relaxation mechanisms. In particular, under the as- 
sumption of an optimal electron spin-spin contact, the 
role of microwave power and electron-nucleus interaction 
was investigated. 

The numerical results presented in Section |III| point 
out how the hypothesis of partial saturation introdocued 
in [19 , not only improves the agreement between TM 
theory an the experimental data of steady state nuclear 
polarization, but also predicts a more realistic behaviour 
for the frequency dependence of the nuclear build up 
time. A similar qualitative agreement is obtained by 
mantaining the full saturation assumption included in the 
original Borghini model and relaxing the constrain of per- 
fect electron- nucleus contact, in presence of a weak, elec- 
tron independent, nuclear spin lattice relaxation term. 
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The two options considered, both fairly good in account- 
ing for the behaviour of the nuclear reservoir, generate 
completely different scenario in respect of the electron 
polarization, as widely discussed in Section |IV| in the 
light also of the experimental observations available in 
literature. 

In order to make the reading of the manuscript more 
fluent, several Appendixes collecting most of the relevant 
mathematics have been added at the end of the main 
text. 



II. MODEL OVERVIEW 

A system made of Nn nuclear spins I (ff, where 1 = 1/2 
for simplicity) resonating all at the same Larmor fre- 
quency LOn, and Nf, electron spins S (f) with mean Lar- 
mor frequency (~ 10^a;„) is considered. The electron 
frequency distribution (ESR line) is supposed to be in- 
homogeneously broadened by single ion anysotropy (i.e. 
spreading of g- factors), the dipolar electron- electron in- 
teraction having no influence on it |23] . In this limit 
the ESR line can be conveniently decomposed in a se- 
quence of narrow individual spin packets of frequency 
uji = ciJe — Ai, width Suj and relative weight fi such that 




Frequency 



FIG. 1: Discretization of the inhomogeneously 
broadened ESR line. Each basic packet centered at 
LJi = ujf, — Ai is characterized by a width 5uj and a 
weight fi. The irradiation frequency = uj^ — Aq and 
the nuclear Larmor frequency aj„ are also indicated. 

The system is assumed to be ruled by five processes 
(see [20] for detailed description): 

1. Microwave irradiation: a single electron transition 
within the microwave (MW) irradiated packet zq. 
The characteristic time constant is named Timw- 

t ^ ; 



2. Spectral Diffusion: a transition involving four elec- 
trons that conserves both energy and total polar- 
ization with a characteristic time constant that we 
named T2e- 

;tt; ^ t;;t 

3. ISS process: a transition consisting in a simulta- 
neous flip-flop of two electron spins, belonging to 
packets separated by w„ , compensated by a nuclear 
spin flip. The characteristic time constant is named 

i t ^ - t i ^ 
Wi cji+cjn n Tiss Wi+cjn n 

4. Electron spin-lattice relaxation: a single electron 
transition stimulated by the contact between the 
electron system and the lattice which leads Pe^i to- 
wards the Boltzman thermal equilibrium polariza- 
tion Pq- The characteristic time constant is Tie. 

t or i ^ t ^ 

I-Pq 
+ 2 
Tie W,; 

5. Nuclear spin-lattice relaxation ('leakage'): a sin- 
gle nucleus transition promoted by the contact be- 
tween the nuclear system and the lattice (with- 
out electron assistance) which leads P„ towards the 
Boltzman thermal equilibrium polarization Po,n ~ 
0. The characteristic time constant Tin is generally 
longer than 10'^ s. 

^ or ^ ^ ^ i 
n Tin n 

A detailed description of the system would also require 
considering transitions between nuclear spins (spin dif- 
fusion) and computing the local dipolar interaction as a 
function of the specific spatial spin location (for both spin 
diffusion and ISS processes). Such approach leads to a 
system of N^ + Nn interacting non linear equations, quite 
difficult to manage. To simplify the problem, we adopt 
a mean field approach where the rate Tiss describes, as 
an effective parameter, both the nuclear spin diffusion 
process and the interaction of two generic electrons be- 
longing to packets i and i + Sup (being 6np the num- 
ber of packets corresponding to aj„) and a generic nu- 
cleus n. Electrons of the same packet are set identical by 
definition and characterized by a local polarization Pe,i, 
whereas a unique polarization P„ and inverse tempera- 
ture Pn is assigned to the whole nuclear system: 

P„(t) =tanh[/3„(i)<5np)]. (1) 
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In this way the number of variables is reduced from Ne + 
Nn to Np + 1. 

The system so defined is studied in the TM regime, 
where the spectral diffusion processes, mediated by the 
electron dipolar interaction, are far more efficient than 
any other process {T2e ^0). In this limit, even when 
the system is out of equilibrium because of the MW ir- 
radiation, a unique spin temperature is established at 
all times among the electron packets (Appendix |A]) . The 
electron polarization Pe.i{t) can thus be written as: 



Pe,.(t) =tanh[/3e A,- c(t))]. 



(2) 



where c(t) and the unique inverse temperature /3e(t) are 
time-dependent parameters. 

The dynamics of Pe,i{t) and consequently of P„(t) is 
determined not only by the highly efficient spectral dif- 
fusion but also by the remaining four processes. Their 
effect is described by the system of rate equations intro- 
duced in and here reported for convenience of the 
reader (under the assumption of finite rates for all the 
four processes). 



dt Tie TlMW 

27iss 

dPn{t) _ Pon - Pn{t) 



(3) 



dt 



Tlr 



where 6i,io is a Kronecker delta and n_ = n_(i, t), 11+ = 
n+(i,i),n„ = n„(i,i) are given by the expressions: 

n- - Pe^^-Snp{t)-PeAi)+Pn{t) [l~Pe,^-Snp{t)Pe,^{i)\ 
n+ = Pe,^+Sn,{t)-P,^,{t)+Pn{t) [l-Fe.»+5n, (t)Pe,« (0] 
n„ = Pe,,+5n,{t)~PeAt)+Pn{t) [l-Pe,,+5n, (t)Pe,« (i)] • 

For numerical computation a discrete time step dt is 
introduced: 



dt 



1 



(4) 



where W^mw = NJq/Timw, We = NjTie. Wiss = 
NeY.hh+SnjTiss and Wn = Nn/Tir.. After each el- 
ementary evolution step according to Eq.(|3|, the ef- 
fect of spectral diffusion (acting on a typical time scale 
5t w ^ dt) is accounted by imposing that the polar- 
izations Pe^iit + St) satisfy Eq.([2]) and the conservation 
of the energy and total polarization: 



n[Pe4t+st)-Pe4t)]^o 



E 



/,A, [Pe4t+st)-Pe4t)] = 0. 



(5) 



In this work we investigate three distinct regimes 
where, in addition to spectral diffusion, one or two more 
processes are assumed infinitely efficient. 



A. Regime I ('Borghini') 

The evolution of the system is derived under the fol- 
lowing assumptions: 

• Tiss — J> 0: a perfect contact between the electron 
and the nuclear reservoirs which allows to estab- 
lish a common electron-nucleus inverse tempera- 
ture /3(i) = (3e{t) — Pn{t) at all times (see Ap- 
pendix |B]) ; 

• TiMW — > 0: a full saturation of the irradiated 
packet io which corresponds to assume Pefi(t) — 0, 
so that c{t) = Aq. 

The steady state solution Pe,i{t -> oo), Pn{t — ?> oo) can 
be computed by solving numerically the know Borghini 
relation (Appendix |C 1[ ) . The dynamics of electron and 
nuclear polarizations can be obtained from the system of 
rate equations ([3|, conveniently adapted to this regime 
(Appendix |D 1| . Moreover, it is possible to write the 
rate equation for the inverse temperature /3{t) (Appendix 
ID 1|), whose solution is not an exponetial function. 



B. Regime II ('partial MW saturation') 

The evolution of the system is derived under the fol- 
lowing assumptions: 

• Tiss — 5> 0: a perfect contact between the electron 
and the nuclear reservoirs which imposes, as in 
regime I, ^{t) = ^^it) = Mt); 

• TiMW 7^ 0: an incomplete saturation of the irradi- 
ated packet io. 

The steady state solution is now function of two variables 
/3 and c and can be evaluated by numerically solving a 
system of two equations, which are a generalized version 
of the Borghini relation (Appendix C2 and [H]). The 
evolution of Pe.i{t) and Pn{t) can be estimated by means 
of the syste m of rate equations ([s]) adapted for this regime 
(Appendix |D2[). 



C. Regime III (poor electron-nucleus contact) 

The evolution of the system is derived under the fol- 
lowing assumptions: 

• TiMW — > 0: a full saturation of the irradiated 
packet «o which imposes c — Aq: 

• Tiss ^ 0: a poor contact between the electron 
and the nuclear reservoirs, modulated by the cor- 
responding parameter l/Tissj which leads (in pres- 
ence of leakage) to two different inverse tempera- 
tures for electrons and nuclei, i.e. Pe{t) 7^ Pn{t)- 
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The steady state solution is now function of two vari- 
ables Pe and Pn and can be evaluated by numerically 
solv ing a syste m co mposed by the two equations |C2| and 
04] (Appendix [CS]). The solution for Pe^^(i) and Pn{t) 
can be obtained from the system of rate equations (|3]), 
conveniently adapted for this regime (Appendix |D 3[ ). 

In the limit Tigs ^ Tie, the contact between the elec- 
trons and the lattice is more efficient than the contact 
between electrons and nuclei. The steady state polar- 
ization profile Pe,i is then achieved in a typical time of 
the order of Tie independentely from any feature of the 
nuclear reservoir. As a consequence, the nuclear system 
'sees', through Tiss, an electron thermal bath at constant 
temperature and the rate equation for P„ [t) assumes the 
linear form: 

P'Jt) = A-BP,,{t) (6) 
where A and B are constant terms defined as: 



B 



On 



21isS^^n ^ J- In 

■ ^ ^ fifi-\-Snp (l Pe,i-\-5np-Pe,i) 



2TissNn 



In 



being Pe^^ given by the Borghini relation (Eq.(C2) in ab- 
sence of nuclei). Its solution: 



Pn{t) = -[l-eM-Bt)] 



(7) 



is an exponential function with a steady state P„ — A/B 
and a polarization time Tpoi = 1/B. 



III. NUMERICAL RESULTS 

The three regimes introduced in Section |ll] have been 
explored by computing a set of build up curves (i.e. po- 
larization versus time) for different values of the signifi- 
cant parameters: zq, PimWj Tiss and Tin- All the other 
parameters of the rate equations have been set as fol- 
lows: Nn/Ne = 1000, Tie = 1 s, A^p = 15, Sup = 3 and 
fi defined according to a Gaussian function with a full 
width at half maximum Awe = 63 MHz and truncated 
at 3(7. This set of parameters is choosen to represent a 
sample of [l-13C]-pyruvic acid doped with 15 mM trytil 
radical in a magnetic field Bq = 3.35 T, at temperature 
T = 1.2 K . Such well known mixture was argued to po- 
larize via TM [TB] and has been studied experimentally 
in great detail [IHl HH [H] . The build up curves obtained 
from the numerical simulation have been fitted by the 
phenomenological law: 



P{t) = Po [1 - cxp i-t/TpoiT 



(8) 



where Pq is the steady state value of the polarization, Tpoi 
is the polarization time and a is a stretching exponent. 
For a = 1, the usual exponential function is recovered. 



A. Regime I: Borghini model 

The physical scenario emerging under the assumptions 
defining the regime I is summarized in FIG. [2] Panel 
A shows the steady state nuclear polarization P„ as a 
function of the microwave frequency wmw in absence of 
leakage and when Ti„ = 10000 s. The same curve is 



obtained by solving Eq.(C2). The calculated values are 



generally higher than those experimentally observed es- 
pecially when moving from the centre to the edges of the 
ESR line. A maximum nuclear polarization of 0.85 is 
reached in absence of leakage when the microwave fre- 
quency is set to ojMW.opt = We - 43 MHz (corresponding 
to the irradiation of the packet io =4). Leakage has only 
a moderate effect on the curve, leading to a 13 % reduc- 
tion of the maximum polarization level for Ti„ = 10000 
s. 

Another interesting quantity for comparison with ex- 
periments (see Section IV) is the average electron polar- 
ization (Pe(i)), defined as: 



(Pe(i)) =^/,;Pe,.(t). 



(9) 



The steady state value (Pg) of this quantity is reported 
in panel B of FIG. [2] When the irradiation frequency 
^MW is close to uje, the ESR line is effectively saturated, 
i.e. (Pe) ss 0. Conversely when the irradiation frequency 
is set at the edges of the ESR fine, (Pe) — >■ Pq because of 
the low weigth fi of the side packets. 

The dynamical evolution of the spin systems can be de- 
rived by means of Eq.(D2| and Eq.(Dl|. The behaviour 



of the nuclear polarization time as function of ujmw is 
shown in FIG.[2j panel C. The larger is the shift between 
We and ll>mWi the shorter is Tpoi or, in other words, the 
steady state is achieved faster when the edges of the ESR 
line are irradiated. 

Finally, on panel D of FIG. [2] the nuclear and electron 
build up curves are displayed for WMW.opt (i-e- with io 
— 4) and no leakage starting from the thermal Boltz- 
man equilibrium condition P„ w and Pf, i = Pq, Vi. As 
introduced in the previous Section, the same inverse tem- 
perature (3{t) (see inset) characterizes both nuclear and 
electron reservoir. The dynamics is non exponential as 
can be noticed by the mismatch between the best fitting 
curve and the simulated data if a stretching exponential 
a = 1 is assumed. The function describing the dynamics 
of /3{t) is computed in Appendix D 1 Despite nuclei and 



electrons share the same temperature over time, since the 
hyperbolic tangent is a non linear function, the nuclear 
and electron polarization build up times are slightly dif- 
ferent, although both in the order of 10'^ s. The initial 
condition f3{t = 0) « originates as follows: 

• when MW are switched on the packet io is imme- 
diately saturated (Pe,o = 0) due to the assumption 
PiMW 0; 

• the fast spectral diffusion (T2c 0) imposes: 
Pe^t = 0) = tanh [/3(t = 0) (A, - Aq)]; 
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FIG. 2: Overview of nuclear and electron polarization in TM-DNP under the assumption defining regime I (T2e = 
s, Tiss = s, TiMW = s) at 3.35 T and 1.2 K. The nuclei/electrons ratio has been set to Nn/Ne — 1000, whereas 
the electron longitudinal relaxation time is assumed to be Tie = 1 s. Panel A: Steady state nuclear polarization P„ 
as a function of the irradiating frequency wmw in absence of leakage (squares) and with Ti„ = 10000 s (circles). 
Panel B: Average electron polarization (Pg) as a function of wmw in absence of leakage (squares) and with Ti„ = 
10000 s (circles). Panel C: Nuclear polarization time Tpoi as a function of wmw in absence of leakage (squares) and 
with Tin = 10000 s (circles). Panel D: Nuclear and average electron polarization build up curves Pn{t) (squares) and 
{Peit)) (circles) at WMw.opt in absence of leakage. The growing curve of the inverse temperature P{t) (triangles) is 
represented in the inset. The non linearity of the differential equations which regulate this regime is reflected in the 
lack of agreement between the calculated trend of Pn{t) and 13 {t) and the exponential best fittings (solid lines). 



• the effective contact between electrons and nuclei 
gives: /3(t = 0) = /?„(t = 0) «0. 



B. Regime II: partial MW saturation 

An overview of the regime characterized by a partial 
saturation of the ESR line is presented in FIG.|3] In panel 
A and B, P„ and (Pe) as function of wmw are shown for 
TiMW = 0, 0.1 and 1 s, i.e. moving from high to low 
MW power. The effect of a partial saturation is twofold: 
on one side a significant reduction of the maximum nu- 



clear and electron polarization values is observed. On 
the other side a clipping of the wings of the steady state 
nuclear polarization curve occurs. This latter effect leads 
to a better agreement with experimental observations, as 
already pointed out in [TS] . 

The behaviour of the polarization time versus wmw 
(panel C) is, interestingly, completely different from 
regime I. As long as the irradiation frequency is close to 
'^ej Tpoi is relatively short, becoming longer and longer 
on moving towards the edges of the ESR line (from 1800 
to 6800 s when Timw = 0.1 s and from 4800 to 15000 s 
when TiMW = 1 s). As expected, longer is T1MW7 less 
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FIG. 3: Overview of nuclear and electron polarization in TM-DNP under the assumption defining regime II {T2e ~ 

s, Tiss = s, TiMw 7^ 0) at 3.35 T and 1.2 K, with N,jNe = 1000, Tie = 1 s and in absence of leakage. Panel A: 
Steady state nuclear polarization P„ as a function of the irradiating frequency ujmw with Timw = (squares), 0.1 s 
(circles) and 1 s (triangles). Panel B: Average electron polarization (Pg) as a function of wmw with Timw = 
(squares), 0.1 s (circles) and 1 s (triangles). Panel C: Nuclear polarization time Tpoi as a function of wmw with 
Timw = (squares), 0.1 s (circles) and 1 s (triangles). Panel D: Nuclear polarization build up curve P„(i) for Timw 
= 0.1 s at WMW.opt- The mismatch between numerical data and exponential best fitting (solid line) points out the 
non linearity of the phenomenon. The inset shows the corresponding growth of the inverse temperature f3(t). 



effective the polarization mechanism is. 



C. Regime III: poor electron-nucleus contact 



The time evolution of P„ represented in panel D, as 
well as the build up curve of /3 (inset) are similar to those 
obtained in regime I, with a non exponential behaviour (a 
rigorous demonstration is not reported in this case) and 
a typical time constant in the order of 10'^ s. The build 
up of the electron polarization instead is somehow more 
complex and characterized by two different time scales. 
The detail of such behaviour are analyzed in Appendix 

E 



The third regime is characterized by a finite contact 
rate between nuclei and electrons. Nuclear and electron 
polarization were computed for Tjss = 0.1 and 1 s both 
in absence of leakage and with Ti„ = 10000 s. The main 
results are shown in FIG. |4j following the same scheme 
of the regimes discussed above. Panel A shows the de- 
pendence of the nuclear polarization on wmw- In absence 
of leakage, as already discussed in [2^, the contact rate 
Tiss does not affect the steady state but only the dy- 
namics and thus the two curves at different Tjss overlap. 
In presence of leakage P„ is reduced. The longer is Tjss 
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FIG. 4: Overview of nuclear and electron polarization in TM-DNP under the assumption defining regime III {T2e = 
s, TiMW = s, Tiss 7^ s) at 3.35 T and 1.2 K, with Nn/N^ = 1000, Tie = 1 s. Panel A: Steady state nuclear 
polarization P„ as a function of the irradiating frequency wmw with Tjgs = 0.1 s (squares) or 1 s (circles) and in 
absence of leakage (empty symbols) or with Ti„ — 10000 s (filled symbols). Panel B: Average electron polarization 
(Pe) as a function of wmw with Tiss = 0.1 s (squares) or 1 s (circles) and in absence of leakage (empty symbols) or 
with Tin = 10000 s (filled symbols). Panel C: Nuclear polarization time Tpoi as a function of ujmw with Tjss = 0.1 s 
(squares) or 1 s (circles) and in absence of leakage (empty symbols) or with Tin = 10000 s (filled symbols). Values 
represented with empty circles are scaled of a factor of j^^. Panel D: Nuclear polarization build up curves P„(t) 
(squares) at WMW.opt, Tjss = 0.1 s and Ti„ = 10000 s. The solid line represents the best fit to an exponential 
function, while the inset displays the corresponding growth of the inverse temperature /3„(t). 



the higher the decrease is. Moreover, the reduction is 
more significant at the edges of the DNP spectrum, a be- 
haviour that becomes clear in the light of panel C, where 
Tpoi is shown to be strongly increased at the wings of the 
spectrum. As a consequence, nuclear relaxation (with 
rate Ti„ — 10000 s) becomes a strong competing mecha- 
nism with respect to the ISS process, forcing P„ towards 
a lower steady state value. 

Panel B highlights a rather interesting feature of 
regime III: the steady state electron polarization is al- 
most unaffected either by Tjss and Ti„. This indicates 
that the nuclear systems, for sufficiently high values of 



Tiss , is only a spectator of the electrons re-arrangement 
under MW irradiation, playing no active roles in the 
evolution of the electron systems towards its equilib- 
rium. Evolution that proceeds through a two-steps pro- 
cess discussed in Appendix |Ej Nuclei have a 'delayed re- 
sponse' characterized by a time constant in the order of 
Nn/{NeTiss) (about lO'' - 10^ s for the set of parameters 
used here) and by an exponential shape as confirmed by 
the good match between the fitting and the simulated 
data in panel D and demonstrated in Section [iT] The 
exponential time course of Pn{t) stems from the linear 
rate equation (|6]), that further remarks the passive role of 
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the nuclear reservoir in the polarization process of elec- 
trons. Correspondingly the nuclear inverse temperaure 
/3„ builds up (inset of panel D) towards a steady state 
value that, in presence of leakage, is substantially differ- 
ent from the end value of Pe{t)- One has, for instance, 
Pn = 1.84 X l{f s vs (3eit) = 3 X 10* s for Tiss = 0.1 s 
and ri„ = 10000 s. 



IV. DISCUSSION AND CONCLUSIONS 

The original description of the TM mechanism pro- 
posed by Borghini and here analyzed in depth in terms 
of electron and nuclear polarization, polarization times 
and dynamics of the inverse spin temperatures has only 
a partial qualitative overlap with the experimental ob- 
servations. The model overestimates the final values of 
Pn especially at the edges of the ESR line, leading to 
an unsatisfactory shape of the DNP spectrum (P„ ver- 
sus wmw)- In fact, even when the MW frequency is 3 
times (T lower than Wg (approx 80 MHz with our choice 
of parameters) and consequently the electron population 
of the corresponding energy levels is very low, the model 
still predicts a very high enhancement of the nuclear po- 
larization which is quite unrealistic and - more important 
- it is not experimentally observed. Conversely, by relax- 
ing either the constraint of a complete MW saturation 
or the constraint of a perfect electron- nucleus contact, 
lower Pn values and sharper DNP spectrum are obtained 
(panels A in FIG. |3) [I]). 

Furthermore, in the Borghini regime, the dependence 
of the efficiency of the polarization transfer on the mi- 
crowave frequency (Tpoi vs wmw, panel C in FIG. p| 
disagrees with DNP experiments. Johannenson et al. 
measured the build up time as a function of the irra- 
diation frequency for a sample of [l-13C]-pyruvic acid 
doped with 18.5 mM of trytil (cfr. Figure 3 of ref. [22]) 
and showed that Tpoi is relatively short as long as wmw 
is set between the two values corresponding to the posi- 
tive and negative maximum of nuclear polarization (DNP 
optimum frequencies) whilst becoming longer and longer 
on moving towards the edges of the ESR line. Remark- 
ably, the correct qualitative behaviour of the polarization 
time is recovered under the assuptions underlying both 
regime II and regime III, as shown in panels B of FIG. [3] 
and FIG. IH 

Up to here, as long as the nuclear parameters only 
(P„ and Tpoi) are considered, regime II and III are both 
in qualitative agreement with the experimental obser- 
vations and nearly superimposable. Actually the two 
regimes are very different, as can be understood from 
panels B and D of FIG. [3] and FIG. H For limited MW 
power and Tigg = 0, the nuclear system and the electron 
one share always the same inverse temperature, generally 
lower than the achieved /3 in case of full saturation. In 
fact, as shown in panel B of FIG. [3j (Pe) tends to the 
frequency-independent equilibrium value Pq when Timw 
increases, as the competition between the MW pumping 



and the electron spin-lattice relaxation unbalances the 
steady state towards the Boltzman equilibrium. Being 
the electron system weakly affected by MW irradiation, 
it is not anymore a forceful source of polarization for nu- 
clei. 

On the other hand, in case of finite electron-nucleus 
contact and Timw =0, the electron system reaches in 
short time a steady state polarization profile, charac- 
terized by a nearly constant inverse temperature /3e. 
In absence of leakage the nuclear system sees only the 
constant-temperature electron reservoir and, with a char- 
acteristic time dependent on the contact ratio 1/Tiss 
reaches a final inverse temperature /?„ = /Sg. In pres- 
ence of leakage the nuclear reservoir is on one side in 
thermal exchange with the electron system at /3e and, on 
the other side, with the lattice at (3^ (x 1/T. The final 
nuclear inverse temperature /3„ is a trade-off value be- 
tween /3e and /3i and, as well as nuclear build up time, 
depends on the two contact parameters Tjss and Tin- 

In order to discriminate which scenario fits better with 
the experimental observations, data about the behaviour 
of electrons must be considered. A valid attempt to char- 
acterize the electron system was made by Ardenkjaer- 
Larsen and collaborators and it is reported in [ISJ [52] . 
By measuring the shift of the ^'^C resonance line (Mi) 
caused mainly by the dipolar fields associated to the po- 
larized paramagnetic centres, the authors indirectly es- 
timated the average electron polarization (Pg) according 
to 111^: 

Ml = ^7re7e7n^^e(Pe) (10) 

where ^ is a coefficient which depends on the shape of 
the sample, 7e is the electron gyromagnetic ratio, 7„ is 
the nuclear gyromagnetic ratio and is the number of 
electrons per unit volume. 

In particular, in [22] the dependence of the nuclear 
shift, and thus indirectly of (Pg), on the MW frequency 
and power was measured. The average electron polar- 
ization (Pg) at a fixed MW power was found to depend 
on ujmw with a behaviour superimposable to that re- 
ported in panel B of FIG. [2] |3j |4] where the degree of 
electron saturation is higher for wmw = Wg — (5a; and 
lower when moving towards the edges of the ESR line 
and, as expected, when the MW power is reduced. It 
is worth to notice that Mi increases rapidly at low mi- 
crowave power and then reaches a plateau for power on 
the order of 40 - 60 mW. For a direct comparison of ex- 
perimental and computational data, the dependence of 
simulated levels of P„ and (Pg) as function of the MW 
power, expressed by Timw, for different values of wmw 
is reported in FIG. [5] The same qualitative behaviour is 
obtained in experimental and calculated data. In numer- 
ical simulations the plateau is reached for Timw i£ 0.05 - 
0.1 s, whereas experimentally a plateau of Mi is reached 
above few tens of mW. Such values are lower than the 
power level commonly used in DNP experiments at low 
temperature (T « 1.2 K), thus suggesting that the as- 
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FIG. 5: Steady state polarization of the electron (panel A) and nuclear (panel B) spin systems as a function of the 
microwave power, expressed by the parameter Timwi for "^mw — "^MW.opt (circles) and wmw = — (squares). 
Remaining parameters are set as follows: Tif, = s, Tjss = 0, Tie = 1 s, Ti„ = cx), Nn/N^ = 1000, Sup — 3, Np = 15. 



sumption of full saturation is more appropriate than the 
hypothesis of partial saturation in interpreting DNP re- 
sults collected in this temperature range. In [19], Janin 
et. al. argued that the increase of the microwave power 
could lead to a heating of the thermal bath which com- 
petes with the polarizing action of the MW themselves, 
affecting the equilibrium nuclear polarization and ending 
up in lower steady state P„ values. Such a heating effect 
however, would also affect the equlibrium electron polar- 
ization, contributing positively to the total saturation of 
the ESR line. In that case, (Pg) should go to zero on 
increasing the MW power (1//3 oo), instead of going 
to the low temperature plateau observed in [22] ■ Thus, 
the assumption of partial MW saturation, although suc- 
cessful in improving the description of DNP from the 
nuclear point of view, shows an intrinsic weakness in ac- 
counting for the electron behaviour. The model of finite 
electron-nucleus contact on the other hand, has similar 
capability in describing the nuclear system, but without 
explicitly contraddicting the experimental behaviour of 
electrons. Overall, given the low temperature DNP ex- 
perimental data available so far, relaxing the condition 
Tiss — appears more promising than removing the sat- 
uration condition Timw = 0. 

An elegant experimental test for better judging the 
physical meaningfulness of regime II and III would 
consist in measuring the electron polarization profiles 
fiPe.ii^Mw)- The expected trends for the two regimes 
are shown in FIG. |6] for ujmw = i^MW.opt (panel A, C) 
and for wmw ^ oje ^ (panel B, D) and described by 
Eq.([2]): under the assumption of regime II /3e = /3„ and 
no packets are fully saturated, whereas in regime III in 
absence of leakage /3e Pn and the irradiated packet is 
characterized by = 0. Especially when ojmw is set 
close to We the electron profile of the two regimes are 
considerably different. 

In summary we have presented the articulated picture 
of thermal mixing DNP generated by the five parame- 



ters model introduced in [20 , in the limit where = 0. 
Three cases in particular have been discussed in detail: 
the Borghini regime, characterized by a strong saturation 
of the ESR line and by a perfect contact between elec- 
trons and nuclei (Timw and Tjss — 0) , the regime of par- 
tial saturation of the ESR fine (Timw and Tjss = 0) 
and the regime of finite electron-nucleus contact (Timw 
= and Tigg ^ 0). The former regime has been shown 
to be less accurate in accounting for the available experi- 
mental observations, whereas the latter two are both ca- 
pable of properly capturing more features of the nuclear 
spin dynamics, whilst predicting different behaviour for 
the electron system. Additional dedicated experiments 
would be desirable in order to clarify which of the two 
predictions gives a better picture of the physical reality, 
although the finite electron-nucleus contact regime looks 
more consistent than partial saturation in describing ex- 
periments at variable MW power. 



There are facts which are not explained within the pic- 
ture exposed in this work, such as the inverse depen- 
dence of the nuclear steady state polarization on elec- 
tron concentration, observed systematically when the ra- 
tio Ne/Nn exceeds a certain value. The statistical ap- 
proach introduced in [20 however, can be extended to 
explore regimes with limited efficiency of the electron- 
electron interaction or, in other words, with limited ther- 
mal contact between different electronic packets. More- 
over, the approach is flexible enough to allow the intro- 
duction of additional interaction terms. By exploiting 
these residual opportunities, we are confident that also 
the still unexplained behaviours will find suitable inter- 
pretation within the general framework of thermal mix- 
ing. 
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FIG. 6: Steady state electron polarization Pe^iit — ;> oo) as function of the electron frequency w in regime of partial 
saturation (Timw = 1 s, solid line) and of non perfect electron-nucleus contact (Tjgs = 1 s, dashed line) for wmw = 
WMW.opt (panel A) and for wmw — uJe — (panel B). Weighted electron polarization fiPe,i{t — >■ c») as function of 
the electron frequency lo in regime of partial saturation (Timw = 1 s, solid line) and of non perfect electron-nucleus 
contact (Tiss = 1 s, dashed line) for wmw = WMW,opt (panel C) and for wmw ^ uje — Suj (panel D). Remaining 
parameters are set as follows: = s, Tie = 1 s, Ti„ = oo, Nn/N^ — 1000, Sup — 3, Np — 15. 
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Appendix A: Spin temperature in the Thermal 
Mixing regime 

Abragam and Goldman [TT| gave a description of TM 
DNF based on the separation between electron Zeeman 
and non-Zeeman contributions in the magnetic Hamilto- 
nian {J^) of the system: 

i i 

with .y/fz accounting for the (Zeeman) interaction of elec- 
trons with the external magnetic field Bp, Jifpfz describ- 
ing the (non Zeeman) interaction of the electrons with 



the local field induced by the anisotropic (^-factors distri- 
bution and representing the nuclear (Zeeman) inter- 
action with Bq. 

Starting from the system Hamiltonian one may de- 
rive the energy of a single electron spin St (belong- 
ing to packet i) associated to its two possible states: 
up (t) of energy eJ — h/2{uje — A^) and down (l), 
Ef = -h/2{uJe - Ai). When MW are off the system 
is at thermal equilibrium with the lattice, at an inverse 
temperature /3l — h/{2kBT) (where fcs is the Boltzman 
constant) and the probability for the spin S'* to be in the 
state up is given by the Boltzman weight: 

p\ oc exp [-/^^(we - Ai)] . 

When MW are on, the system is out of equilibrium. If 
now the existence of a unique temperature among the 
different packets is postulated as in [11], the probability 
pj can be expressed in terms of a generalized Boltzman 
weight: 
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oc exp 



nA,, 



2kBTa ^ksTp 



= cxp [-{auje - /3Ai)] 



where the two parameters a = h/ {iksTa) and /? = 
h/ {2kBTi3) are normally referred as Zeeman and non Zee- 
man inverse temperature respectively. The polarization 
of the spin can be then written as: 



P, = 



Pi - Pi 
Pi +Pi 



— — tanh [a We — /3 A^] 



The same espression can be derived by observing that, 
whenever a process much faster than the other events 
ruling the system exists, the detailed balance for such 
a process must be satisfied at any point in time. In all 
the TM scenario considered in this work, the 'spectral 
diffusion' mechanism depicted here below 

has been always assumed to be a fast process. Its cor- 
responding detailed balance condition can be written in 
terms of the fraction of electrons up - P^^ [t) - and of the 
fraction of electrons down - P~^{t) - at time t : 

Pt-8{t)PeAt) ^ PU*)Pe,,+sit) 

Then, by using the relation 



PtM = 



l + PeAt) 



PeAt) = 



l-PeM 



2 ' ' 2 

one comes to an equation for the electron polarization 

{I - Pe,i){l + Pe,^-S) ^ {I + Pe,j){l - Peg+s) 
(1 + Pe,«)(l - Pe,^-S) (1 - Pe.j){l + Pe.j+s) 

which is satisfied if 

Pe,.(t) =tanh[/?e(i)(A^-c(t))], 

qed. 



an equation for the nuclear polarization is obtained: 

Pe.iit) — Pe,i+Snp{t) 



Pn{t) 



1 ~ Pe,i+Sn„{t) Pe,i{t) 



(Bl) 



Being Pea{t) = tanh [^(t)(Aj - c(t))], Eq.(|BT]) can be 
rewritten as: 

Pn{t) =tanh [l3{t)Snp], 

which defines the existence of a unique common tem- 
perature between nuclear spin system and electron non 
Zeeman reservoir. 



Appendix C: Steady state equations 

Following jTT], the energy of the whole electron and 
nuclear systems can be split into a Zeeman electron con- 
tribution 

and into a non-Zeeman electron plus Zeeman nuclear 
term 

ENZ-zAt) = ^N„huJnPn{t) - ^ A^e^ ^ A,/,Pe,^ W ■ 

The time evolution of the two energy reservoirs is de- 
scribed by the following equations: 



dt 2 



E_f ^0 ~ Pe,i{t) r Pe,0 
J I Tf. Jo 



Tie 



dt 2 



^A,/,%^-f/oAo 



TiMW 
Pe,0 



Tie 



TlMW 



2 " " Ti„ 



(CI) 



Since at the steady state the two reservoirs have con- 
stant energy by definition, simultaneous vanishing of 



both right-hand sides of Eq.( CI ) can be imposed to eval 



uate the steady state polarizations Pe,i and P„. 



Appendix B: Detailed balance for the ISS process 



1. Regime I: Borghini 



When the ISS process is assumed to be as fast as spec- 
tral diffusion (regime I and II in this work), also the de- 
tailed balance for such a process has to be verified at all 
times. Thus, the following condition needs to be satisfied: 

it) P;,^+8n, it) Pn it) = P-^ it) P++Sn, it) P^ (0 • 

Using the relations 



p:^it) 



l + Pnit) 



Pnit) 



l-Pnit) 



For regime I, being c = Ag, a unique equation in- 
dipendent on Timw is required to evaluate /?. Thus, by 



multiplying the first Eq.(Cl ) by Ap/oje and adding it to 



the right-hand side of the second Eq.(Cl), one gets rid 
of the microwave transition probability I/Timw and ob- 
tains a single expression, which is the celebrated Borghini 
relation: 



MA, - Ao)Pe.. + AqPo - UJn^P^Pn = 0. (C2) 
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2. Regime II: Partial saturation 

For regime II, both the two equations are required to 
define the unique solution for /3 and c: 



Pq — Pe,i 
Tie 



fo 



Pe,0 
TlMW 
Pe,0 







-tie JlMW J Ira 



0(C3) 



3. Regime III: poor electron- nucleus contact 



The Borghini relation [C2] holds also in this regime, but 
it is not sufficient to determine unambiguously /3e and /3„. 
The additional relation which completes the steady state 
picture is the rate equation for Pn{t) that, once imposing 
the stationary condition, writes: 



JVe 



2TissiV, 



(^1 Pe,iPe,i-\-6np^ ~\~ 



(C4) 



Appendix D: Electron and nuclear spin dynamics 

The numerical procedure described in the main text 
to evaluate Pe,i{t) and Pn{t), when all the processes but 
the spectral diffusion have a finite transition rate, has 
been conveniently adapted for the three regimes consid- 
ered. The strategy consists in using conservation laws to 
manage all mechanisms assumed to be infinitely efficient, 
while computing rate equations only for processes with 
finite rate. 



1. Regime I: 'Borghini' (Tiss and Timw — > 0) 

Rate equations are used to account only for the effect 
of the electron and nuclear spin-lattice relaxation: 



dt Tie 
dPnit) Pon~Pn{t) 



dt 



Tin 



(Dl) 



time of the transitions T2e, Tiss and Timw- These equa- 
tions are conveniently written as: 

- 51 [PeAt+St)-PeAt)] = /O [Pe,0{t+St)-Pe,0{t) 

^ /,A,; [PeAt+StyPeM] +/0A0 [PeAt+St)-Pe,0{t) 
-SPMW^ -^OJn [Pn{t+St)-PAt)] = 

SO that, by means of simple algebra, the condition: 

^/,(A,-Ao) [PeAt+St)-PeAt)] 



Fast processes (spectral diffusion, electron-nucleus con- 
tact and MW saturation) are accounted by the conserva- 
tion of the total polarization (for variations induced by 
spectral diffusion or ISS) and of the total electron non 
Zeeman plus nuclear Zeeman energies: 

^ /. [PeAt + 5t) - PeAt) - duoSP^'"^] = 

^ /,A, [PeAt+6t)-PeAt)] -^^n [PAt+6t)-Pnit)] =0 

where JP^^^ indicates the variation due to MW irradi- 
ation and the time step i5t — being the characteristic 



N, 



[PAt+st)-PAt)]^o 



(D2) 



is obtained. By solving this equation one derives P{t+dt) 
and computes Pe,i{t + St) — ~ tanh [f3{t + 6t){Ai — Aq)] 
and PAt + St) = tanh [/3(i + St)Snp]. 

It is interesting to study also the evolution of the 
inverse temperature f3(t) that in regime I, as demon- 
strated in the Appendix |B] is the same for both the 
electron non Zeeman and the nuclear Zeeman reservoirs: 
f3{t) = f3e{t) = Pn{t)- Moreover, since full saturation 
imposes c = Aq, I3(t) is the only unknown variable of 
the problem. Hence, by means of Eq.(D2) and Eq.(Dll, 



it is possible to describe analitically the time behaviour 
of /3{t). At a generic time t + dt, by assuming St — > 0, 
Eq.(|D2l) writes: 



h {- tanh [I3{t -h dt) (A, - Aq)] - P,,,(t + dt)) 
N 

(A,-Ao) -^^n [tanh (/3(i + dt)uA - PA^ + dt)] = 



-Po - tanh [/?(<) (A, - Ao)] | (A, - Aq) 
-ujn [tanh {f3{t + dt)ujn) — tanh (/3(i)(u;„)] 



Now, using Eq.(Dl) for replacing Pe,i{t + dt) and P„(i- 
dt) one obtains: 

h {- tanh [/?(< + dt) (A, - Ao)] + 
-ftanh [I3{t) (A, -Ao)] + 
dt 

_Nn 
= 

Then, with the first order expansions: 
I3{t+dt) w f3{t)+p'{t)dt 

tanh [/3(t-fdi)x] wtanh [l3{t)x] +P'{t)x {l-tanh^ [f3it)x]} 

and some algebric calculations, the following equation for 
/3{t) is achieved: 

13' {t) {J2 - ^o)' [1- tanh^ [/3{t) (A, - Aq)]] + 



-ujI [l-tanh^ [(5{t)iOn]] 



(D3) 



Ti 



^ ^ /, (A, - Ao) {Po + tanh [/3(i) (A, - Aq)]} = 
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Eq.(|D3|) is conveniently rewritten as: 



E fi (Ai - Ao) tanh (A, - Ao)] - AqPo 



iVeTi, 



-w2 [l-tanh' im^n]] 



m 

after negleting with good approximation the term 
Y: h (A^-Aof [1- tanh^ [m (A,-Ao)]] . 



Although no temptative of solving analically Eq.(D4) 



is done here, it is clear that the solution can not be an 
exponential function, as already anticipated in the main 
text. 



2. Regime II: partial MW saturation {T2e, Tiss 0) 

The system of rate equations is used to describe the 
effect of partial MW saturation as well as electron and 
nuclear spin-lattice relaxation: 



P, 



e,0 



IMW 



dt Tie 
dPnit) ^ Pon-Pnjt) 
dt Tin 



whereas spectral diffusion and electron-nucleus interac- 
tion are accounted by the following conservation laws: 

^/,A, [p,At+styPeAt)] 

N 

-^^n [PAt+st)-PAt)]^o 

with dt being the characteristic time of the tran- 
sitions and Tiss- By solving this system one ob- 
tains /3{t + 5t) and c(t -I- 6t) and computes Pe,i{t + 
St) = -tanh[P{t + St){A,-c{t))] and PAt + St) = 
tanh [(3{t + St)Snp]. 



3. Regime III: poor electron-nucleus contact {T2e 

and TiMW 0) 

The system of rate equations takes into account the 
effect of the electron-nucleus contact and of the electron 
and nuclear spin-lattice relaxation: 

dPeAt) ^ Po-PeAt) 
dt Tie 

27iss 

dPAt) Pon - PAt) Ne 



dt 



Tin 2TigsA^„ 



The effect of the other processes (spectral diffusion and 
full MW saturation) is accounted by the conservation of 
the total polarization (when the variation is induced by 



spectral diffusion) and of the total electron non Zeeman 
plus nuclear Zeeman energies: 

h [PeAt + - PeAt) - '5,,,„<5P^^] = 
^/,A, [PeAt+5t)-PeAt)]=^ 

where SP^^^ indicates the variation due to MW irradi- 
ation and the time step St ^ Q being the characteristic 
time of the transitions T2e and Timw- These equations 
are conveniently written as: 

- 51 [PeAt+5i)-PeAt)\ = /O [Pe,o{t+5t)~P,^o{t) 

/.A,; [PeAt+St)-PeAt)] +/0A0 [Fe,o(i+'5i)--Pe.,:(0) 
_SpMW^ = 

so that, after simple algebric calculations, the following 
condition is obtained: 

Y hi\-Ao) [PeAt+St)-PeAt)] = 0, 

that allows deriving f3{t+St) and thus computing Pe.i{t + 
St) = - tanh [(3 {t + St) {A^- Ao)]. 



Appendix E: Dynamical behaviour of the electron 
average polarization 

The evolution of the average electron polarization (Pg) 
in regime II and regime III shows a peculiar behaviour 
characterized by two different time scales, as sketched in 
FIG.E 



1. Partial MW saturation 

In this regime, due to the hypothesis T2e = and 
Tiss = 0, one has PeA*) = tanh [/3e{t) {A, - c{t))] and 
Pe{t) = Pnit) = I3(t). The evolution of {P^) is deter- 
mined both by I3{t) (panel D, FIG. [s]) and c{t). 

At short times {t « Tie) the inverse temperature /3 is 
determined by the large nuclear system for which /3„ (t — 
^) — Pl ~ 0. When Timw — 0, the only solution is 
c{t = 0) = Ao and consequently PeA^ = 0) = 0,Vi. For 
partial saturation (Timw > 0) the profile of Pe.i becomes 
a flat function (corresponding to the condition c{t = 0) — > 
00) which quickly evolves with a characteristic time Tie 
towards an intermediate level between and Tqj that can 



be calculated using the first equation of system C3 



Pe,i — Pq 



T 



foTi, 



IMW 



IMW 



At longer times (i « Tpoi) the evolution of (Pe) is 
mainly due to f3(t) dynamics (being c{t) approximately 
constant) and it is thus characterized by a time constant 
in the order of 10'^ s. 
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FIG. 7: Panel A: Build up curve of the average electron polarization {Pe{t)) in regime of partial saturation (regime 
II) for TiMW = s (Borghini limit, thick solid line), 0.1 s (small-dashed line) and 1 s (large-dashed line). Remaining 

parameters are set as follows: Tie = 1 s, Ti„ = oo, N^/N^ — 1000, Iq =4, 6np = 3, Np = 15. Panel B: Build up 
curve of the average electron polarization {Pe{t)) in regime of poor electron-nucleus contact for Tjgs = s (Borghini 
limit, thick solid line), 0.01 s (dotted line), 0.1 s (small-dashed line) and 1 s (large-dashed line). Remaining 
parameters are set as follows: Tie = 1 s, Ti„ = 10000 s, Nn/N^ = 1000, io = 5, Sup = 3, Np = 15. 



2. Poor electron-nucleus contact 

The dynamics of both f3e{t) and {Pe{t)) is character- 
ized by two time scales: a first rapid component with a 
characteristic time in the order of Tie and a second slow 
component with a characteristic time in the order of Tpoi. 

In this regime, due to the hypothesis — and 
TiMW = 0, one has Pe,i{t) = tanh [/3e(t)(Ai — Aq)], with 
/3e(i) 7^ fin{t)- Depending on the time scale considered, 
the system can be qualitatively depicted and /3e esti- 
mated accordingly. 

• At very short times (t — >■ 0), being the contact be- 
tween electrons and nuclei finite, the electron sys- 
tem is unaffected by the presence of the nuclear 
reservoir and reaches immediately the inverse tem- 
perature Pb predicted by Borghini and defined by 
EqjD2] after setting = 0. 

• After this initial 'thermalization' phase, at times 
t « Tie, the electron reservoir is on one side in con- 
tact with a thermal bath at temperature 1//3b (de- 
termined by interaction with the lattice, by spectral 



diffusion and by the highly effective MWs), while 
feeling on the other side the nuclear ensamble hav- 
ing an initial temperature /3„ = /3l ~ 0. Thus, on 
a time scale of few Tie, the inverse temperature /3e 
moves towards a target value between Pb and /3l, 
depending on the strenght of the two contact times 
Tie and Tjss- When the electron- nucleus contact is 
poorly efficient /3e — /3b; conversely /3e — 5^ /^l ~ 
for strong electron- nucleus contact. 

• At large times [t Tpoi), /3n(i) evolves from (i^ 
towards its final steady state /?„ and fie{t) evolves 
as well, reaching an intermediate value between /?„ 
and j3B- 

In summary, as long as the electron-nucleus contact is 
poorly efficient, the electron inverse temperature /3e is 
only slightly affected by the nuclear reservoir and it is 
thus seen by this latter as a constant value equal lo (3b- 
As discussed in Section |ll] and in Section IIVI this be- 
haviour leads streightforward to an exponential build up 
curve for nuclear polarization. 
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